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Abstract. Differential models for hydrodynamic, passive-scalar and wave turbulence 
given by nonlinear first- and second-order evolution equations for the energy spectrum 
in the fc-space were analysed. Both types of models predict formation an anomalous 
transient power-law spectra. The second-order models were analysed in terms of 
self-similar solutions of the second kind, and a phenomenological formula for the 
anomalous spectrum exponent was constructed using numerics for a broad range 
of parameters covering all known physical examples. The first-order models were 
examined analytically, including finding an analytical prediction for the anomalous 
exponent of the transient spectrum and description of formation of the Kolmogorov- 
type spectrum as a reflection wave from the dissipative scale back into the inertial 
range. The latter behaviour was linked to pre-shock/shock singularities similar to the 
ones arising in the Burgers equation. Existence of the transient anomalous scaling and 
the reflection-wave scenario are argued to be a robust feature common to the finite- 
capacity turbulence systems. The anomalous exponent is independent of the initial 
conditions but varies for for different models of the same physical system. 
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1. Introduction 

Differential models of turbulence, such as the Leith model ([T]) and the Kovasznay model 
([2]) have played an important role for achieving a qualitative and even quantitative 
understanding of turbulence in various physical situations. Differential models have 
been also used in wave turbulence. In particular, for the water gravity wave turbulence 
such a model was introduced in |3j. The main advantages of the differential models 
are their relative simplicity and their great flexibility when it comes to the inclusion 
of new physical effects, e.g. simultaneous presence of a strong and a weak turbulence 
components in superfluid turbulence ([3]), modelling atmospheric turbulence with two 
scale-separated sources ( 0 ), new types of forcing or dissipation like reconnections of 
superfluid vortices ( 0 ) ( 0 ) or sound radiation by vortices ( 0 )- 

Another advantage of the differential models is that they allow to study non¬ 
stationary turbulence. Let us consider an isotropic hydrodynamic turbulence with an 
initial spectrum in a finite range of wave numbers k = |k|. How will such a system 
evolve toward the Kolmogorov power-law spectrum with a ^-independent flux of energy 
in the inertial range of scales? One could guess that the spectrum would spread out 
of the initial finite range of scales toward high values of k leaving the Kolmogorov 
spectrum immediately behind the moving front. However, recent findings obtained 
within the integro-differential and differential closures of turbulence, as well as the 
exact results for the Burgers equation, indicate that this may be an incorrect view. 
By solving numerically the kinetic equation for weak MHD turbulence, it was found in 
[8] that the initial spectrum behind the moving front is indeed a power law, E ~ k~ x , 
but it has an exponent ~ 7/3 which is clearly distinguishable and steeper than the 
one of the Kolmogorov-type spectrum (2 for the weak MHD turbulence). It was 
found that the Kolmogorov-type spectrum gets formed in the system by a reflection 
wave propagating from high to low wave numbers after the spectral front reaches the 
dissipation scale. Similar behaviour was observed for the differential (Leith) and integro- 
differential (EDQNM) closures of Navier-Stokes turbulence in [9] and in respectively. 
Rigorous bounds for the anomalous index in the Leith model were recently found in 
HU: the index was proven to be strictly greater that Kolmogorov’s 5/3 and strictly less 
than 1.95 which agrees with the numerical value ~ 1.85. 

For the Burgers equation, it was noted long ago in [12j that the typical wave 
breaking scenario in which a pre-shock (cubic root singularity) is followed by a shock 
(jump discontinuity) in the k -space corresponds to an initial formation of a power-law 
spectrum with exponent 8/3 which is gradually replaced by the Burgers x = 2 spectrum 
invading the fc-space as a wave moving from high to low k. 

Presently, it is understood that the described above scenario in which an anomalous 
power-law precursor spectrum appears prior to the Kolmogorov-type spectrum must 
be a generic property of the finite capacity systems, i.e. the systems for which 
the Kolmogorov-type spectrum is integrable at k = oo (meaning that for vanishing 
dissipation the Kolmogorov-type spectrum contains only a finite energy). This scenario 
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is described, using the classification of ra, by a self-similar solution of the second kind. 


The fact that a turbulent system belongs to the finite capacity type does not depend on 
the turbulent closure and, therefore, we should expect that the exact system exhibits 
the same anomalous scaling behaviour without any assumptions or approximations 
related to a particular closure. However, the levels of resolution and precision of the 
current experimental measurements and numerical simulations are currently insufficient 
to identify a scaling range prior to formation of the Kolmogorov spectrum which would 
be wide enough to be able to distinguish confidently between the Kolmogorov and 
the anomalous exponents. On the other hand, the anomalous power-law index, being 
independent of the initial conditions, does depend on the closure used, as will be seen 
in the MHD and Navier-Stokes turbulence examples later in this paper. Whether there 
is a universal value for the anomalous index and if such a value could be traced to 
particular singularities of the dynamical fluid equations (like in the Burgers equation 
example) remain questions to be answered. 

In the present paper we extend the study of the anomalous power-law to a general 
three-parametric class of the differential models, which is relevant to several known 
examples of strong and weak finite capacity turbulence systems. In addition to the 
second-order equations of the nonlinear diffusion type generalising the model of [T], we 
study simpler first-order nonlinear wave equations generalising the turbulence model of 
[2]. Those second-order models can be solved exactly. They yield an anomalous power- 
law scenario, that may be explained by the formation of Burgers-type singularities - 
pre-shocks and shocks. 

The paper is organised as follows. In section [2J we describe a generalized version of 
Leith Model that can be used as a toy model for turbulent systems with a finite-capacity 
spectrum. In section |3j the calculations of [9] and [IT] are extended to infer the existence 
of self-similar solutions whose energy spectra are steeper than Kolmogorov’s. In section 
[4j we use numerics to characterize this anomaly, and construct a phenomenological 
fitting formula for the anomalous exponent as a function of the parameters of the 
model. In section [5j we discuss an even simpler class of models which include the 
Kovasznay model of 3D turbulence [2j, and for which anomalous exponents exist and 
can be analytically derived. Finally, in section [7] we present a summary and discussion 
of our results. 

2. Nonlinear diffusion models for turbulence 

We consider the following inhomogeneous non-linear diffusion procfirst-orderess in k- 
space for the energy spectrum E(k,t), 



( 1 ) 


where d = 1, 2 or 3 is the dimension of the system and indices m and n are real (usually 
rational positive) numbers depending on the physical system. Constant C in general 
has a physical dimension which depends on the particular system. 
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The model ([!]) is constructed so that it has two fundamental power laws E ~ k~ x 
as stationary solutions - a constant flux Kolmogorov-like spectrum with x = Xk = 
(m — d)/{n + 1) and a thermodynamic equilibrium spectrum with x — xt — 1 — d. 

One can check that the diffusion model also admits the general time-independent 
solution 

E s {k) = k d ~ l [Ak 1 - rn+{ ' 1 - d)n + B) 1/[1+n) (2) 

with arbitrary positive real constants A and B. Case B = 0 corresponds to the 
Kolmogorov-type spectrum and case A = 0 to the thermodynamic spectrum. Thus, 


the case with both A and B having non-zero values is a mixed state which was named 
a “warm cascade” in [9]. Note that at k —> 0, the mixed state ([2]) is dominated by 
the Kolmogorov (flux) part and for k —> oo, this solution tends the thermodynamic 
equilibrium spectrum. 

Note that equation (Jlj) describes a conservative model, J 0 °° E dk = const. Naturally, 
viscosity can be taken into account by adding to the right-hand side of this equation 
term —vk 2 E, where v = const > 0 is a viscosity coefficient. Other types of dissipations 
may also be modelled depending on the physical system. However, for the second- 
order models we will only be concerned with non-st at ionary initial states of evolution 
at which the dissipative scales are not reached and, therefore, the dissipative terms can 
be neglected. (We will be able to have a more general treatment including the stage 
where the viscosity is important when we consider the first-order models). 

Also, in this paper we will be concerned with finite capacity systems for which the 
Kolmogorov-type spectrum is integrable at k —> oo. In the other words, no matter how 
large the dissipation wave number is, the Kolmogorov state will have a finite energy 
density in the physical space. Thus, we are interested in the physical examples with 
xr = (m — d)/(n + 1) > 1. Let us now mention some examples. 

2.1. Examples 

2.1.1. 3D hydrodynamic turbulence The case with d = 3 ,m = 11/2 and n — 1/2 
corresponds to the Leith model introduced in |T]. In this case constant C is dimensionless 
and has value ~ 1/8 ([§]). 

2.1.2. Passive scalar turbulence In this case E has the meaning of the passive scalar 
spectrum and d — 2 or 3. Also, since the advection is passive, the resulting equation 
must be linear in E, therefore n — 0. The constant m depends on the roughness of the 
advecting velocity field. For smooth velocity (Batchelor’s regime) we have m = d + 1 so 
that Xk = 1 (|14j) . Actually, the Batchelor regime is on the borderline of the infinite 
capacity systems and strictly speaking our analysis will be inapplicable to this case. 
However, we will see that the case xk = 1 + e is interesting and it will nicely illustrate 
a smooth transition between the finite and the infinite capacity regime asymptotics. 

For a rough velocity corresponding to the Kolmogorov spectrum (Obukhov-Corsin 
regime) we have m = d + 5/3 so that xk = 5/3 ([15| 16]). 
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2.1.3. Wave turbulence In this case d — 1,2 or 3 and n is the order of the resonant 
wave interaction less two, e.g. n = 1 for three-wave processes, n = 2 for four-wave 
processes, etc. Indeed, the nonlinearity degree of the differential equation is n + 1 i.e. 
the same as the one of the respective integro-differential wave-kinetic equation. The 
constant m depends on the particular type of the waves. Particular examples include: 

(i) Gravity water waves: In this case d = 2, n = 2 (a four-wave processes) and 
m = 19/2 so that Xk = 5/2 (see |17j : for the differential model for this system 
see [18] and [IS]). 

(ii) Capillary water waves: In this case d — 2, n — 1 (a three-wave processes) and 
m = 11/2 so that x K = 7/4 (see e.g. [2U]). 

(iii) Sound waves: In this case d — 3, n — 1 (a three-wave processes) and m = 6 so that 
x k = 3/2 (see [2.1]; the differential approximation for this system is the same as 
the one for isotropic MHD turbulence derived in [22]). 

(iv) Alfven waves: In this case d = 2 (because even in 3D the system is very anisotropic 
and the spectrum is almost two-dimensional) n — 1 (a three-wave processes) and 
m — 6 so that Xk — 2 (see [8]; for the differential model for this system see m 

(v) Kelvin waves on vortex filaments: In this case d — 1, n — 2 (a four-wave processes) 
and m — 6 so that xk = 5/3 (see [24] and [25]). 

3. Self-similar solutions of the nonlinear diffusion models 

We are interested in the early evolution of the spectrum which is initially non-zero only 
in a finite range of k. As the front of the spectrum propagates toward large k, for the 
values of k which are much greater than the initial /c’s the solution tends to a self-similar 
solution of the second kind. It is to be found by making in equation ([Tj) a substitution 
in the form 


E(k, t) — (£* — t) a F(r]), r] = k/k in and fc* = c(t* - t)A (3) 

Such a solution describes an explosive propagation of the spectral front k = fc* in a finite 
time t*, i.e. fc* —» oo as £* —» oo (function F(rj) is zero at r) > 1). Thus, the solution 
exists only for a finite time after which one can no longer ignore viscosity to describe 
the subsequent evolution. 

Obviously the constant c depends on the initial spectrum: the stronger turbulence 
is the faster it evolves. Solution for the case with c / 1 can be obtained from the 
solution for the case with c = 1 by a rescaling. Thus, thereafter we will put c— 1. We 
will also put C — lm equation ([!]) because it can be absorbed into the time variable. 
To make the self-similar framework consistent, or in other words to make the equation 
for F involve the similarity variable ij only, a and (3 need to be chosen so that 


na + (m — d — l)/3 + 1 = 0 . 


( 4 ) 
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For a steady state power-law spectrum, E ~ k~ x the self-similarity implies F ~ p~ x . 
The meaning of such power-law asymptotics will become clear shortly. For now, we 
express a and f3 in terms of x = —a//3 using the condition Q: 

« =- , i - and I 3 = - , i -• (5) 

nx + a + 1 — m nx + d + 1 — m 


We then obtain the following equation for the profile function F: 


1 

d + 1 — m + nx 


[xF + pd v F\ = d v \rf n F n d ri (F/'q d ~ 1 )] . 


( 6 ) 


The special case corresponding to the Leith model, d = 3,m = 11/2 and n = 1/2, was 
treated in !)] and in nu • In this case, l/(d + 1 — m + nx) = 2/(x — 3). Below, we will 
extend the analysis to the general finite capacity case. 

Following [9] and HU, we realise that any initial spectrum concentrated in a finite 
range of k will asymptotically tend to a self-similar solution for large k. This corresponds 
to a solution of Equation ([6]) with two boundary conditions: F(l) = 0 corresponding to 
a sharp propagating front of the spectrum and F ~ p~ x at rj <C 1 corresponding to a 
stationary power-law spectrum forming behind the propagating front. This formulation 
is a nonlinear eigenvalue problem because it has a solution only for one value of x = x*. 
It is precisely the “eigenvalue” x* and its dependence on m, n and d which we will aim 
to find below. 

Like in [9] and HU, we will first find an autonomous system equivalent to equation 
(J6|. In addition to being autonomous, the new system must have non-singular fixed 
points. This is achieved by choosing parametrically a suitable new time variable t(F, rf) 
. The latter can be found first using the following change of variables, 


F = p^ v - x and F'{rj) = A rj-^p^a, 


(7) 


where p and A must be selected so that one gets a suitable autonomous system (see 
Appendix [B|, namely : 

p — l/n and \ = (m — d—l)/n. (8) 


It is now convenient to work with the new independent variable r = t(t}), which we 
define according to the formula 


dr 1 
d log rj p ’ 


so that 


d 

d7? 


1 d 
pp dr 


(9) 


Plugging these changes of variables into Equation 
autonomous system : 


© then yields the following 


p\ T ) = -p(p + 


a r = 




d — 1 

p 


+ d - 2^j (M) 
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Pi = (0,0), P 2 = 0, 


h 


A (a: — A) 


and P 3 = (1, -1) 




(A — 1)(A + d — 1) 


( 12 ) 


The nature of the fixed points can be established by examining the Jacobian matrix 
A(p, er) of system (10) - & see appendix [Aj We find that P\ is a saddle-node and 


P 2 is a saddle. The nature of P 3 depends on the parameter values. Namely, P 3 is an 
unstable node for x < x_, an unstable focus for x G [x_; x c ] a stable focus for x G [x c ; x+] 
and a stable node for x > x+, where x± and x c are explicitly defined in formulae ( |A.4 ) 
and (A.3). A local Hopf bifurcation of creation of a limiting cycle around P 3 occurs 
at x — x c . As found in Oj and ra the vicinity of the point Pi corresponds to the 
rj — > 0 part of the solution whereas point P 2 corresponds to the sharp front, rj — 1. The 
goal is to find such x = x* that one could have an orbit connecting P\ and P 2 , i.e. a 
heteroclinic orbit. Such a heteroclinic orbit arises at a global heteroclinic bifurcation 
at x = x+ associated with creation of a heteroclinic cycle consisting of two heteroclinic 
orbits one of which being the required solution and the other one is the piece of the 
cr-axis connecting P\ and P 2 (see details in ra¬ 
in the special case of the Leith model, m showed that the value of x* lies above 
xk and below x c . This result still holds true for the general class of models studied 
in the present paper. While this bounding can be used to obtain specific asymptotic 
behaviours for the critical exponent x* (see Appendix [C]) , it does not indicate how this 
exponent precisely varies for “intermediate regimes” of the parameters (m,n,d). Still, 
this theoretical bounding of x* can be used as a first proxy to numerically determine 
x* with a dichotomic search algorithm. Our aim is twofold : i) observe the general 
dependence of re* with the parameters used in the model and ii) provide a practical 
“engineering” approximation for x*(m, n) for the physically motivated values of the 
parameters (m, n, d ) previously described. The algorithm that we use and the numerical 
results are described in the next section. 


4. Numerics 

J h l. Dichotomic algorithm 

The value x*, such that the trajectory starts at P\ and ends at P 2 can easily be 
determined numerically using a dichotomic algorithm. The basic idea is simple: if 
we start near P 2 and compute the trajectory backwards in r for x x*, we will miss P\ 
and will end up either spiralling into the focus P 3 or crossing the cr-axis to the quadrant 
of positive p’s. Detecting these events will allow us to iterate x closer to x*. 

In the sequel, the set of equations (10) - ([TTj) is rescaled so that P 3 = (1, —1) in 
the (p, cr)-plane. 


(i) Set xi = X- and x 2 — x + . 

(ii) x i (xi T x 2 )/2. 
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Figure 1. Illustration of the dichotomic algorithm for m = 11/2, n = 3/2, d = 3. 
Here e = 10 -2 and ?? = 10 -6 


(iii) Compute Pi, P 2 , P 3 . Choose an initial point Pj in the vicinity of P 2 , say 

a = -p 2 + 


(iv) Integrate equations (10) - (11) backwards in time until either cr(—r) > 0 or 
|P(—t) — P 3 I < e (e <C 1 is prescribed). This sets a final time 77 . 

(v) If cr(-Tf) > 0 , then x 2 <5— x else X\ x. 

(vi) Repeat steps (b) to (e) until |x 2 — Xi| < 77 (7 is a prescribed accuracy of the result). 

(vii) Return x* = (xi + x 2 )/2. 

For m = 11/2 and n = 3/2, the algorithm yields x* ~ 1.85; see Figure [1} 


4-2. General Behaviour. 

Recall that the value of (xk, n, d) determines the triplet (m, n, d) as m — d — (n + l)x/<. 
The dichotomic algorithm was run, in order to observe how the value of x* varies with 
xk , w and d. For each set of parameters, the value of x* was determined with a precision 
of rj = 10 ~ 5 . 

We considered 50 different values of n’s logarithmically spaced between e ~ 3 ~ 0.05 and 
3, and 50 values of xk — l’s logarithmically spaced between e ~ 5 and 6 . With this choice 
of parameters, we aim to investigate the behaviours n —» 0 + which describes the passive 
scalar and Xk —> 1 + , the limit below which the energy spectrum has no longer a finite 
capacity. Note that we found that the dimension d had no significative quantitative 
influence on x*, and therefore only discuss in details the numerical results obtained for 
the case d = 3. 
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A 3D rendering of the general behavior of the deviation of the critical exponent x* 
from the Kolmogorov value xk is displayed in Figure [2j The deviation is here measured 
by the quantity Sx = (x* — xk)/xk > 0 and plotted against xk and n. The figure 
indicates the following trends. 

(i) For fixed n, \xk — x*| —> 0 as xk —> 1. 

(ii) For fixed n, \xk — x*| —* oo as xk —> oo. 

(iii) For fixed xk, \xr — a?*| —>• 0 as n —* oo. 

(iv) For fixed xk, \xr — x*| —> oo as n —> 0. 

Note that the asymptotic behaviors (a) and (c) observed in the numerics can be 
proven mathematically using the bounding theorem of mi. It suffices to recall that 
X- < Xk < x+ < x +J explicitly compute x + — x_, and observe that this difference 
vanishes in both asymptotics. On the contrary, the behaviors (b) and (d) cannot be 
obtained with the same argument, as in the latter case the difference x + — X- diverge - 
see Appendix [Cj 

It is interesting that the limits Xk —> 1 and n —> 0 do not commute. Recall that 
the case n — 0 corresponds to the passive scalar turbulence and, in particular, the 
case (n = 0 ,xk = 1) corresponds to the passive scalar turbulence in Batchelor regime 
(smooth velocity). For the latter regime it is actually known that Xk = x * in the 
case when turbulence is forced at large scales, i.e. the Kolmogorov-type (Batchclor- 
Kraichnan) spectrum develops right behind the propagating front; see |26j . In absence 
of forcing the spectrum is non-universal. This is not surprising as this case describes a 
(border-line) infinite capacity system. However, it is rather striking that in the passive 
scalar turbulence with rough velocity field x* —> oo. This can be interpreted as a 
spectrum that decays faster than any power law. In other words, the Kolmogorov 
regime (Corrsin-Obukhov) gets established as a backward wave propagating from the 
dissipative wave numbers to the smaller wave numbers. There is no precursory scaling 
in the larger wave numbers in this case, i.e. in the opposite direction with respect to 
the infinite capacity systems’ behavior. 

4.3. An engineering fit for the physical range of parameters. 

A closer analysis of our numerics indicate the existence of n-dependent power-law 
behaviors for x* as a function of the Kolmogorov exponent xk in both asymptotics 
Xk —> 1 + and Xk 1, namely : 

x* — 1 ~ eP°( n \x K ~ 1)^^ as log(xA- — 1) —>■ ±oo. (13) 

For the values of n here considered, these power-law regimes seem to hold for xk > 
xfi = 1.4 and xk < xfi = 1.1 (see Figure [3). Note that the behaviors of the /3^i(n)’s 
(Figure |3| is compatible with the asymptotics described in the previous section. Indeed, 
for large n for example, we observe that fif ~ 1 and (3$ ~ 0, implying x* ~ Xk a feature 
indeed in agreement with the limit (n —> oo) previously described. 
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Figure 2. 3D rendering of the deviation from the Kolmogorov exponent, measured 
through the quantity 5 = (a;* — xr)/xr plotted against Xk and n for the case d = 3. 
The same plots for the cases d = 1 and d = 2 would lie slighlty below, but would be 
practically undistinguishable from the case d = 3. 


In a more practical perspective, we also wish to get a ready-to-use proxy for the 


anomalous exponent x* associated to the physical examples described in Subsection |2.1 
The behavior [13] suggests to model the logarithm of x* — 1 as a n-dependent piecewise 
linear function of log(xiy — 1). We therefore propose the following “engineering fit” : 


log (x* - 1) = < 


Pi ( n ) log (xjf - 1) + P 0 ( n ) 

Pi (n) log (x K ~ 1) + Po(n) 
If-l 


if x k < x K = 1.1 
if xk > x\ = 1.4 


K + 7 T- jf (log ( x k ~ 1) - Ir) otherwise, 


/+ _ / 
"K l K 

± 


(14) 


with l K = log(x x - 1), /*=/?! log(x x - 1) + A 




Po in) = exp a oi lo g * n , 


and j3f I 


n) = 


1 + exp 'Y af { log* n. 


i =0 


i=0 


The a^’s are fitting coefficients whose values are given in Table [l] These numbers 
reproduce the numerical data with less than 0.5 percent of error on the value of x* 


for the range of parameters here considered. However, the formula (14) should not 


be extrapolated boldly outside its range of validity, in particular in the regions of the 
parameter space where it is expected that the deviations from Kolmogorov regime are 
very large - namely here Xr > 7 or n < 1.05. The 4 th order polynomial fittings of 
the /Ts indeed give a n —> 0 behavior which disagrees with the asymptotics previously 
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logiofZif - 1) 



Figure 3. Left : Illustration of the piecewise-linear engineering formula (14). The dots 
are the values obtained with the dichotomic algorithm, while the lines correspond to 
the formula |l4| . For representation purpose, not all the data points are displayed. The 
color codes n (see right panel). The blue dotted lines separate the regions Xk < x^, 
Xk < Xk < %k an d x k > x f- Indication of a n-dependent linear behavior is apparent 
for x < x]( and x > xf. Right : Evolution of [3 1 (n) and /3 0 (n) as measured from the 
dataset (see Formula (14) for the definitions). 


described. 


5. The first-order model 


One of the drawback of the generalized Leith models (|Tj) is that we do not know whether 
the anomalous spectrum can be entirely and explicitly determined mathematically. Let 
us now consider an even simpler model of turbulence based on a nonlinear first-order 
PDE: 


d t E = —Cdk(k p E q ) — vk 2 E, 


(15) 


where real (usually positive rational) constants p and q are chosen so that in the case 
of zero viscosity coefficient v there is a power-law constant-flux solution of Kolmogorov 
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d=3 


d=2 



i=0 


i=l 


i=2 


i=l 

1 

i= 

4 


a 0i 

-1.14- 

10 1 

-1.35- 

10 1 

-6.14 

o 
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-1.405 

o 

o 

-1.31 

10“ 

-1 

a u 

-1.28 ■ 

10 1 

-1.34- 

10 1 

-6.05 

o 

o 

-1.36 

o 

o 

-1.27 

io- 

-1 

a 0i 

-2.79- 

o 

o 

-2.92 

o 

o 

—1.15 

o 

o 

-2.98 ■ 

o 

1 

-3.48 

10- 

-2 

a H 

-2.83- 

o 

o 

-1.78- 

o 

o 

-1.30 

o 

o 

-3.79 • 

o 

1 

-4.93 

10- 

-2 

a 0i 

-8.46 • 

o 

o 

-9.98- 

o 

o 

-4.26 

o 

o 

-9.02 • 

o 

1 

-7.78 

10- 

-2 

a u 

-9.88- 

o 

o 

-9.91 

o 

o 

-4.20 

o 

o 

-8.61 • 

o 

1 

-7.34 

10- 

-2 

a 0i 

-2.49 ■ 

o 

o 

-2.77- 

o 

o 

-1.21 

o 

o 

-3.46 ■ 

o 

1 

-4.22 

10- 

-2 

a ti 

-2.85- 

o 

o 

-1.35- 

o 

o 

-1.16 

o 

o 

-3.59 • 

o 

1 

-4.87 

10- 

-2 


Table 1. Fitting constants used for the formula (14) in the cases d = 3 and d = 2 


type E ~ k~ XK , i.e. p = q xx■ Note that this model no longer has a thermodynamic 
equilibrium solution. The second relation between p and q and the (usually dimensional) 
constant C depends on the particular physical problem. 

Let us mention again our main examples and specify the values p and q for the 
particular systems. 


5.1. Examples 

5.1.1. 3D hydrodynamic turbulence The case with p = 5/2 and q = 3/2 corresponds 
to the Kovasznay model introduced in [2]. In this case constant C is dimensionless. 

5.1.2. Passive scalar turbulence In this case E has the meaning of the passive scalar 
spectrum and, since the advection is passive, the resulting equation must be linear in 
E; therefore q — 1. The constant p depends on the roughness of the advecting velocity 
field. For Batchelor’s smooth-velocity regime p = 1 so that Xk = i am)- For a rough 
velocity corresponding to the Kolmogorov spectrum (Obukhov-Corrsin regime) we have 
p = 5/3 so that xk = 5/3 (||I5] H6j). 

5.1.3. Wave turbulence In this case q is the order of the resonant wave interaction 
less one, e.g. q — 2 for three-wave processes, q — 3 for four-wave processes, etc. The 
constant p depends on the particular type of the waves. Particular examples include: 

(i) Gravity water waves: In this case q — 3 (a four-wave processes) and p = 15/2 so 
that xk = 5/2. 

(ii) Capillary water waves: In this case q — 2 (a three-wave processes) and p = 7/2 so 
that xk = 7/4. 
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(iii) Sound waves: In this case q — 2 (a three-wave processes) and p = 3 so that 
x k = 3/2. 

(iv) Alfven waves: In this case q = 2 (a three-wave processes) and p = 4 so that xk = 2. 

(v) Kelvin waves on vortex filaments: In this case q = 3 (a four-wave processes) and 
p = 5 so that xk = 5/3. 


5.2. Analysis of the first-order model 


Let us make in equation (15) the change of variables E = e x ^ q k p / q (e is the energy 
flux). We have: 

-d t e = —Ck v l q e x ~ x l q d k E — uk 2 £. (16) 

q 


5.2.1. Hydrodynamic and wave turbulence (Case q > 1). 
q > 1. Making a further change of variables 


First, let us consider the case 


U = £ l ~ l ' q 


and £ = k 1 ~ p/q , 


we have 


d t u = C di (w 2 ) — v(q — 1 )i 2 ^ 1 p ^ q) u, 


(17) 


(18) 


where C = C(p — q)/ 2. In all our examples with g / 1, we have 1 < q < p < 3q. 
We later assume that this relation holds. This implies C > 0, and xk > 1 so that 
p/q — 1 > 0 in the aforementioned change of variables ( [17] ) (recall that Xk = p/q). The 
upper bound on p will be useful to describe the dynamics after the shock. 

It is easy to see that in the case v = 0, we recover an inviscid ID Burgers equation. 
Thus, at the initial evolution stage when the viscosity effect is negligible we should 
expect features of the Burgers behaviour, in particular wave breaking leading to a pre¬ 
shock singularity. Thereafter, we can also anticipate a shock formation. However, the 


viscous term in equation (18) is obviously different from the one of the Burgers equation 


and, therefore, the shock structure should be carefully re-examined. 

Starting with the inviscid stage, we first note that (unlike the case of the second- 
order models) the initial conditions for the spectrum E with finite support in k are not 
appropriate: if the velocity at the boundaries of the support is zero then the support 
will not grow in time and there will be no cascade to larger fc’s. Thus we will consider 
the initial spectrum which extends to infinite k (although decreasing there as will be 
discussed later). For simplicity, let there be a minimal wave number k m i n below which 
the initial spectrum is zero. General smooth profiles like this have both positive and 
negative derivative parts, and their evolution will lead to wave breaking with an infinite 
derivative, ie dpa = oo, at a single point £ — 7* at a time t = t *. Since C > 0, the 
breaking will happen somewhere where dqu > 0. We will suppose that the initial profile 
U()(£) is such that in the range of positive 8(Uq{£) the second derivative is negative, 
dfju 0 (£) < 0. In this case the breaking occurs at Z* = 0. Provided that <9^i o (0) f 0 (this 
corresponds to E(k) ~ a t large k initially), we have u{£) ~ £ 1 ^ at time 
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t — £* in the vicinity of the breaking point for the profile. In terms of the spectrum, 
such a behaviour implies 

E[k) = u^ =l k~ v ^ q ~ £3( q -D~^ — k~ x *. (19) 


where 


x* 


P-Q P 

39(9-1) 9’ 


( 20 ) 


Note that this spectrum is steeper than the Kolmogorov-type spectrum, x* > Xk = p/q- 
Again, we see that the limits p —> q and q —> 1 do not commute. Recall that the 
limiting limit, p = q corresponds to the passive scalar turbulence in smooth velocity 
field—a (marginally) infinite capacity system which is considered later in this section. 
To obtain, the large k behavior of the spectrum, we take into account uq{£) = (k p / q E) q ~ l 
and di = The condition duu 0 (£) < 0 then reads: 


d k (k? /q d k (k? /q E) q ~ 1 ) < 0. 


Condition diUo(0) 0 and the condition that the flux is zero at k = 00 leads to the 
asymptotic behaviour uq(£) ~ £ at £ —> 0. This implies a power-law spectrum E ~ k~ x 
at k —> 00 with x = — r. 

q-l 

For example, the anomalous transient exponent is ./;* = 2.1111 for 3D 
hydrodynamics and x* = 7/3 for the MHD values. Note that the latter value coincides 
with the value obtained by numerical simulations of the kinetic equation in jH]- In both 
cases, we obtain x = 3. 

Now let us consider dynamics after the spectral front reaches the dissipative scale, 
t > £*. This will be characterised by a shock in the profile u(£) near t = 0 the structure 


of which is determined by the dissipation term in equation (18). Within the shock one 
can neglect the time derivative term : 

2Cd l u-u(q-l)( 2 ^ l ~ p / q l = 0 . 


Solving this equation and taking into account the condition u(0) = 0 yields : 


39 -p via — li 

u{£) = Ci - C 2 £ i-p , with C 2 = ~^ 7 Z -r > 0 

C(3q - p) 


(recall that p < 3q is here assumed). C\ > 1 is independent of £ but may be dependent on 
t. The first term on the right-hand side here is negative and singular at £ = 0. Therefore, 
the dissipation makes the spectrum have a cut-off wave number k v corresponding to some 
£ = £ v such that u(£ v ) = 0, i.e. 


C 1 C(3g-p) 
% - 1) 


q-p 
3 q-p 


and k v 


C\C{f3q — p) 
v{q - 1 ) 


q 

3 q~P 


( 21 ) 


The value of C\ is to be determined from the matching to the jump in u{£) arising from 
the inviscid Burgers solution. In the vicinity of the wave breaking point, the velocity 
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profile behaves as £ + u(t — t*) ~ u 3 , so that Cj ~ \Jt — t* for the jump at £ = 0. For 
the spectrum, this implies : 

E(k) = (Ci - C 2 k 3 ~ p/q ) ^ jfc"?. (22) 

i _ E 

For k < k u the spectrum is of Kolmogorov type and satisfies : E(k) = C( k * ~ 

1 _ p 

(t — t*) 2 ('?- 1 ) k i We therefore see that the Kolmogorov-type spectrum invades the k- 
space propagating from large to low Us. It gradually replaces the anomalous spectrum 
E ~ k~ x * (whose amplitude is almost time-independent). 

5.2.2. Passive scalar turbulence (Case q = 1) Let us now consider the case q — 1, 
p > 1, which corresponds to a passive scalar in a rough velocity field. The case can be 
treated by the same method and similar results are found. As before, there exist no 
power-law asymptotics is formed before the dissipative scale is reached. The equation 
for the energy flux now reads : 


d t £ = —Ck p dk£ — vk 2 £. (23) 

Using the further change of variables t = k l ~ p . it becomes 

d t £ — {p — l)Cdk£ — vC^pe. 


Writing e(k, t ) = F(k, t ) with k = £ + (p — 1 )Ct, we obtain 


d t F = —v(k — (p — 1 )Ct) 1 -pF, hence In F = 


v 


3 -P 


-(k- (p-i)cty-p + g(k), 


( P — 3)C 

with G being an arbitrary function. Plugging the original variables back in yields : 


E(k, t) = k-”e~ irao*’ P G(k 1 ~ p + (p - 1 )Ct), 


(24) 


with G{k) = e G ^ - a function defined by the initial condition. For the finite capacity 
case (p > 1), the initial spectrum propagates toward low values of £ = /c 1_p , therefore 
towards high Us. The propagation speed is constant in the space of variable t. This 
means that the front reaches l — 0 (or alternatively k = oo) in a finite time. Therefore, 
for any viscosity u, no matter how low, the dissipative scale 


K 



l 

3 -p 


(25) 


is reached in a finite time f* ~ k(y p /(p — 1 )C. 

For t < U, the evolution is inviscid and non-universal: the spectral slope has 
Kolmogorov value — p near the maximum of function G(k), but it gradually varies as 
one moves away from this maximum. At t ~ £*, the maximum of G(k) reaches k v 
and stays there thereafter, gradually becoming flatter on its left side, k < k v . This 
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corresponds to spreading of the Kolmogorov scaling E ~ k~ p from the dissipative scale 
k u into the inviscid range k < k u —a picture we have anticipated. 

Let us now consider the case q = p = 1, which corresponds to the example of the 
passive scalar turbulence in a smooth velocity field. The energy flux satisfies 

d t e = —Ckd k e — uk 2 £. (26) 


Let £(k,t ) = F(a,t) with a(k,t) = Ck exp(—Ct). We obtain : 


d t logF = —z/ — exp(2Ct), so that F(a, t ) = G(a) exp 2 C 3 
C 2 


(27) 


The function G(a) has to be fixed by the initial condition E(k, 0) = Eo(k). This gives: 


E(k, t ) = E 0 (ke 


-Ct\ 


exp 


uk 2 


(28) 


From this solution we see that the inviscid stage lasts until the time ~ 1/C when the 
front reaches the dissipative wave number k v = \j2Djv. At this stage the spectrum is 
simply being stretched exponentially in time to larger fc’s in a self-similar way conserving 
the total energy. For t > 1/C, this stretching continues in the range k < k u getting 
more and more flat, whereas at k ~ k v the spectrum experiences an effective cut-off. 
Note that at no point one observes the Kolmogorov-type scaling, which is a natural 
behaviour for the free decay in infinite capacity systems. 


6. Numerics 


In order to illustrate the previous analysis, we now take the example of Alfven wave 
turbulence, and show the outcomes of numerical simulations of both modelings: first 
and second-first-orderorder. In principle, the first-order model is fully integrable by the 
method of characteristics. However, the resulting solutions are typically implicit and 
numerical solutions remain helpful to visualize the evolution. 

The details of the numerical simulations are the following. For the second-order 
model, we solve Equation [l] with d — 2, n — 1 and m — 6 and a small viscosity 
v = 5 • 10~ 5 as performed in [2B). The initial spectrum is E(k, 0) oc /c 3 e -fc2//25 , such that 
Eft = 0) = 1. The dynamics is then evolved until time tf — 0.02 using a Crank-Nicolson 
scheme with adaptative time-stepping. 200 values of ’Afs are chosen, logarithmically 
spaced from k — 1 to k — 2 20 . For the first-order model, we use Equation (18) with 
q = 2,p = 4 and add a very small viscosity, namely v — 5 • 10~ 4 for a resolution of 
4,000 points. As an initial condition, we chose Uq(£) = sin(£/4) for 0 < £ < 47T and 
uo(£) = 0 for £ > 47t. This initial profile satisfies duUo(£) < 0 in the range of positive 
d(U Q (£) together with <9yuo(0) f 0. It therefore matches the conditions of the previous 
analysis. Crank-Nicolson and Adam-Bashforth numerical schemes are used to reach the 
final time tf — 2000, choosing dt = 10~ 2 for the time stepping. Note that the initial 
spectrum at large k behaves as E(k) ~ feW -1 )/!? -1 ) — ^-3 w hich is steeper than both the 
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k 



t/u 


Figure 4. Left: Compensated MHD spectra obtained numerically from the second- 
order differential model at different times: earlier than t* (in blue) and later than f* 
(in red). The spectra are normalized by the total energy. The times are reported by 
the dots on the right inset. The blue dots correspond to t/f* — 1 ~ —6 ■ 10 _1 , —4 • 10 _1 , 
—10 _1 , —4 • 10 —2 , —10 —2 , —2 • 10~ 3 , 0; the red dots to i/i* — 1 ~ 4 • 10~ 3 , 4 • 10 -2 , 
— ICC 1 , 4 ■ 10 -1 , 7 ■ 10~ 2 . The right inset shows the time evolution of the energy, 
i* ~ 8.3 • 10 -3 marks the onset of energy dissipation. The figure is made using the 
data obtained by E.Buchlin and previously reported in |23] (in a different form). 


Kolmogorov-type spectrum and the transient anomalous spectrum (see Figure [7]). The 
results for the second-order model are shown on Figure [4j For this example, recall that 
the anomalous exponent is x* ~ 2.088 while the Kolmogorov exponent is Xk = 2. The 
difference between the anomalous transient exponent and the Kolmogorov exponent is 
therefore very small, less than 5 percent. However, this difference is clearly apparent in 
our numerics. Following the evolution of the total energy, one can identify £* ~ 8.3-1CT 3 
as the onset of energy dissipation (see the right panel of Figure |4]). For times earlier 
than i*, the figure clearly displays an anomalous growth of the energy spectrum as k~ x * 
with x* ~ 2.088 as predicted by our analysis. After the spectrum reaches the dissipative 
scale, it is gradually replaced by the Kolmogorov-type spectrum E(k) ~ k~ 2 (so that 
the compensated spectrum behaves as : k x *E(k) ~ fc 0 ' 088 ). This numerical experiment 
provides a clear illustration of the calculations presented in section [3j 

The results for the first-order model are shown on Figures ms and [7} The time 
evolution of the profile u(£, t ) is shown on Figure[5j The figure displays the salient feature 
of the first-order model, namely the formation of a sharp cut-off at small i through a 
wave breaking process. The compensated plots of In figures [6] hint at the dominance 
of the K 1 / 3 inertial range scaling at the times close to t*, and 1° inertial range scaling 
at later times. In terms of the spectrum, this should translate into the development of 
a transient spectrum with index x* = 7/3, gradually replaced with a Kolmogorov-type 
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t 



Figure 5. Left panel : function u(£,t) for the first-order differential model of MHD 
turbulence at different moments of time. Right panel: zoom at the small-£ region. The 
initial profile is in dashed black. The times before the shock are t = 80, 400 (in blue). 
The times after the shock are t = 640, 800 (in red). Please see also the inset on the 
right of Figure [7] 




i £ 


Figure 6. Compensated velocity profiles for the first-order differential model at 
different times. The spectra are compensated by l 1 / 3 (left panel) and 1° (right panel). 
The times displayed are those of Figures[5] The insets show the same data using log-log 
coordinates. 


spectrum with index x = 2 at later times. Such a trend is detectable on the spectra 
shown on Figure [7J corresponding to different times chosen both before and after the 
front has developed. 

7. Conclusions 

In the present paper, we have analysed differential models for various examples of 
hydrodynamic, passive-scalar and wave turbulence given by the second-order Leith-type 
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k 


Figure 7. Spectrum of MHD turbulence obtained numerically from the first-order 
differential model at different moments of times. The times and the color legend are 
the same as in Figure [5] . The right inset shows the time evolution of the total energy, 
with the blue and red dots indicating the times picked in Figures [5] to [7] 


nonlinear diffusion and the first-order Kovasznay-type nonlinear transport equations. 
Both types of models predict formation an anomalous (steeper than Kolmogorov) 
transient power-law spectra. 

The second-order models were analysed in terms of self-similar solutions of the 
second kind, and the anomalous exponent was found numerically and presented via a 
phenomenological “engineering fit” formula for a broad range of parameters relevant for 
the physical examples, both the known ones and potentially the new ones. A peculiar 
non-commuting limit was noted for the systems close to the example of a passive scalar 
in a smooth velocity field. 

The first-order models were examined analytically, including finding an analytical 
prediction for the anomalous exponent of the transient spectrum and description of 
formation of the Kolmogorov-type spectrum as a reflection wave from the dissipative 
scale back into the inertial range. Even though the first-order models are less realistic 
than the second-order models, their advantage is that they admit a full analytical 
treatment, which allows to relate the anomalous transient scaling and the subsequent 
Kolmogorov-type spectrum to the pre-shock/shock singularities similar to the ones 
arising in the Burgers equation. On the other hand, somewhat counter-intuitively, the 
first-order model is much harder than the second-order model to deal with numerically. 
This is due to the appearance of shocks, as typical for hyperbolic systems. 

It is tempting to think that in more realistic settings, including the second-order and 
the integro-differential models (e.g. EDQNM, kinetic equations), a similar link between 
the anomalous transient scaling and a singularity could exist. Moreover, such a link may 
exist also on the dynamical level beyond the closures. With a bit of imagination one 
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Model 

Kolmogorov 

Leith 

Kovasznay 

3D hydrodynamics 

5/3 

~ 1.851 

19/9 ~ 2.111 

Passive Scalar (smooth) 

1 

undefined 

undefined 

Passive Scalar (rough) 

5/3 

oo 

oo 

Gravity Waves 

5/2 

~ 2.514 

11/4 = 2.75 

Capillary Waves 

7/4 

~ 1.799 

2 

Sound Waves 

3/2 

~ 1.509 

5/3 ~ 1.667 

Alfven Waves 

2 

~ 2.088 

7/3 ~ 2.333 

Kelvin Waves 

5/3 

~ 1.672 

16/9 ~ 1.778 


Table 2. Anomalous exponents obtained for various turbulent systems described in 
terms of the (generalized) Leith and Kovasznay described in the paper. 


could suggest that the long-hunted for Euler singularity, if exists, may be responsible 
for the anomalous scaling which is a precursor to the Kolmogorov spectrum. Clearly, at 
this moment in time such suggestion is a pnre speculation. However, we would like to 
recall numerical simulations of [27] (and more recently [28]) where a precursor power- 
law with a steep exponent ~ 4 was found for the hydrodynamic turbulence, and it 
was linked to transient pancake-like quasi-singular structures, ft would be interesting 
to study this effect numerically further using modern high-resolution codes, and also 
check the scenario of formation of the Kolmogorov spectrum as a backscatter wave in 
the Navier-Stokes turbulence. 

Our final remarks are about the degree of universality of the transient anomalous 
exponents. On one hand, these exponents are universal in the sense that they are 
insensitive to the initial spectrum. However, clearly the exponents are model-dependent, 
e.g. they are different for the first-order and the second-order models of the same physical 
systems. For example, for the Navier-Stokes turbulence x* ~ 2.111 for the first-order 
model, x* ~ 1.851 for the second-order model, x* ~ 2 for the EDQNM closure, and 
x* ~ 4 in DNS. Thus, the robust feature common to the different models is the fact 
that the anomalous scaling exist, but not the predicted value of its exponent. 
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The stability of the fixed points Pi, P 2 and P 3 defined by Equation (12) can be 
determined by computing the Jacobian matrix A(p, cr) associated to the system (10) 
-CD- The latter reads 


/ 


— (2p + a) 


p,x 


\X(x — A) 


+ 2(d — 1) 


'1 1' 


y J— 
A p 

P + 


d- 1 
A 


+ d- 2 


P 


a 


A 


x — A 


— 2Acr+ 


~P 

s 

d- 1 
P 


\ 


+ d- 2 


P 


In particular, at the fixed points I \, P 2 and P 3 , the Jacobian matrix respectively takes 
the values 


A x = 


A 3 — 


P 


x — A 

/ 


^ 0 0 ^ 

x/X 1 


/ 


Ao — 


P 


x — X 


0 

x + d — 2 d — 1 

+ —- -1 


\ 


V a 


AyU 


fix 


+ 


X 

~l 

p 


-p3 


and 

A 

~P3 


\X(x — A) 
with p 3 defined as p 3 = 




P3 


P 


x — A 


+ 


P 


(A — 1)(A + d — 1) ’ 


P 

2A — 1 + (d — 1) ( 1 + — 

P 


(A.l) 


For all known examples p — 1/n > 0. Then P 3 is a saddle-node and P 2 is a saddle. To 
determine the nature of P 3 , one needs to determine the eigenvalues of A 3 . We compute 
the trace T 3 and the determinant D 3 of A 3 as: 


T 3 


^ + ( A ( 2 + ^) + (<i_1) ( 1 + ^)" 1 )'’ 3 and D3 ~ P3 ' 


(A.2) 


The determinant D 3 is positive. Indeed, from the definitions of Xk on one hand and 
of (A, pi) on the other hand (Equation ([8])), we obtain that A = (p + 1 )% — fi, with 
fi = 1/n > 0 and xk > 1 for the finite capacity scenarios that we here consider. As the 
dimension d is greater than 1, we conclude that both A, A + d — 1 and hence D 3 are 
positive. The trace T 3 is positive for x < x c and negative otherwise, with x c defined by 
the formula 


x, 


— X- 


(A - 1)(A + d - 1) 


A(2 + 1 / fi) + (d — 1) (1 + 1 / fi) 
We now compute the two solutions x± of the equation Tf = 


- 1 

4A 3 . We obtain 


(A.3) 


x± — X — 


_ P _ 

p 3 (A(2 + l//i) + (d - 1)(1 + 1/p) — 1) =F 2 ^p 3 


(A.4) 


Therefore, P 3 is an unstable node for x < x_, an unstable focus for x G [x_;x c ] a stable 
focus for x G [x c ]x + ] and a stable node for x > x + . As an example, for the Leith model 
d = 3,m = 11/2 and n = 1/2, we obtain ~ 1.01, x + ~ 2.28 and x c ~ 1.95 (and 
x K = 5/3 ~ 1.67). 
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It is straightforward yet slightly tedious to obtain the system 0)-® from the equation 
© for the profile function. For thoroughness, the missing steps are provided below. 

Equation (10) is a simple consequence of the definition of the variables p,a and r 
as provided by the equations ([7]) and ([9]) (p and a are not independent). Indeed, from 
the definition ([7]) and the chain rule, we obtain 


F'(v) = PP M VO)- 


dr 


-V 


—A—1 


— Ap M p 


pt —X— 1 


d log rj 

= p M_1 p _A_1 (pp(t) — A p) on the one hand, 

= p^V^Ao- on the other hand. 


(B.l) 


For non-vanishing p, equating the latter two equalities yield Equation (10). Deriving 
Equation (11) is more involving. First, one gets from Definition (J9]) that: 

dp p'(r) 


dp 


PV 


dcr d (r) 

and — =- 

dp prj 


(B.2) 


One then needs to plug the definitions (Jt|) into Equation Q for the profile function. 
The l.li.s of (j6j) then becomes 


1 


[xp + Aa] p M 1 rj A . 


d + 1 — m + nx 

As for the r.h.s, it now reads in terms of A, p, n and m as follows, 

r/ m-d-(n+ l)A-yn+l)^-2 j( m - d) (1 - d) p 2 + A 2 ^71 + 1 - d 2 

+A ^ i^rn — d -^ + (1 — d)(n + 1)^ pa + A a'(r) 


(B.3) 


(B.4) 


In order to obtain an autonomous equation for a', one needs to equate the powers of p 
and p that appear on both sides. This prescribes the following relations between (A, p) 
and (m, n), 

pn = 1 and m — d — n\ — 1 = 0, (B.5) 

i.e. relations d8l). With such a choice of A and p, the l.li.s and r.h.s of the profile equation 


- as given by the formulae (B.3) and (B.4) — can be considerably simplified. They 
now respectively read: 


Lhs = /.(v> + V) md 

x — A 

r.h.s = (1 — d) + 1^ p 2 + A 2 d 2 + A ^2 — d + -—— 
Equating both sides and dividing by A yields Equation 0- 


(B-6) 


pa + Ad'(r). 
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We compute Ax = x + — x_ in several asymptotics. When the quantity goes to zero, 
then x* —» xk- We use the notations of appendix [A] to compute 


Ax 


4/i 

D \' 2 {A 2 D 3 


—, where 

4) 

and D 3 


A — A(2 + l//i) + (d — 1)(1 + 1 j if — 1 

_ If _ 

' (A — 1)(A + d — 1) ’ 


(C.l) 


Recall also that A = r(xk — 1) + xk- 


(i) For xk —> 1 (fixed n), we obtain A — > 1, A — > d(l + l//x), D 3 ~ /i/(d(A — 1)). Hence 
Ax ~ 4/i 3/2 (A - 1) 3/2 / (d 1/2 (/r + l) 2 ) -A 0. 

(ii) For xk —> oo (hxed n), we obtain Ax ~ 4fi l 11 ^ 2 {/i + 1)xa'/ ((2 + 1 / p) 2 p — 4) —>■ oo. 

(iii) For n —y oo (hxed xk), we obtain Ax rs j 4p?I 2 [xk — 1 ) 3 / 2 /(xk + d — l) 1 / 2 — > 0. 

(iv) For n —> 0 (hxed xk), we obtain Ax rs./ p 1/2 {x K - 1) -A oo. 

For the cases (a) and (c), we therefore obtain that the deviations goes to 0, viz., 

x* — xk —> 0. For the two other cases, we cannot conclude. All that we can say is that 

the deviations cannot grow faster than xk (case (b)) or ttT 1 / 2 (case (d)). 
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